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Abstract. We present the time integrated and time resolved spectral analysis of a sample of bright bursts selected with F pea k > 
20 phot cm -2 sec -1 from the BATSE archive. We fitted four different spectral models to the pulse time integrated and time 
resolved spectra. We compare the low energy slope of the fitted spectra with the prediction of the synchrotron theory [predicting 
photon spectra softer than N(E) oc E~ 2 ^ 3 ], and test, through direct spectral fitting, the synchrotron shock model. We point 
out that differences in the parameters distribution can be ascribed to the different spectral shape of the models employed and 
that in most cases the spectrum can be described by a smoothly curved function. The synchrotron shock model does not give 
satisfactory fits to the time averaged and time resolved spectra. Finally, we derive that the synchrotron low energy limit is 
violated in a considerable number of spectra both during the rise and decay phase around the peak. 
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1. Introduction 

The nature and emission mechanisms responsible for the 
prompt emission of Gamma-ray bursts (GRB) are still a mat- 
ter of debate. From the phenomenological perspective much 
effort has been made in order to characterize and identify typi- 
cal spectral properties of bursts, mostly by applying parametric 
but very general and simple spectral models. Pro bably the most 
widely adopted is that suggested by Band et al. ( 1993 ), namely 
a smoothly connected double power law model. Its attractive 
feature is that it characterizes, within the observational energy 
window, the most relevant quantities, namely the peak energy, 
representing the energy at which most of the emission occurs, 
and the low and high energy components which are related, 
according to the most accredited emission theories, to the par- 
ticle energy distribution and/or to the physical parameters of 
the emitting region. Indeed GRB spectra result to be well rep- 
resented by the Band parameterization (N(E) oc E a '@), with 
typical low energy power law photon spectral indices a be- 
tween -1.25 and -0.25, high energy spectral indices j3 around 
-2.25 (Preece et al. 1998) and peak energy E p typically around 
1 00-300 keV. 

A fundamental issue in the spectral analysis of GRBs is the 
integration timescale of the spectra which are observed to vary 
even on millisecond timescales (Fishman et al. 1995b). BATSE 



spectra, which are integrated for a minimum of 128ms, there- 
fore represent the best way, presently available, to constrain the 
burst emission mechanism. Furthermore, in order to compare 



the average spectral properties of different bursts, also time in- 
tegrated spectra, covering the duration of the pulse or the entire 
burst, have been used in literature. 

It has been found in general that the Band model not only 
well describes the time integrated spectra, but it also appears to 
fit the time resolved spectra of bright bursts. Alternative spec- 
tral models - a nd in p articular the predict ions for synchrotron 
emission (Katz 1994, Crider at al. 1997a ,b) - have also been 



recently tested on a large sample of t ime re solved spectra by 
Preece and collaborators (Preece et al. 2000). 
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Within this scenario, here we present the study of the spec- 
tral properties of single pulses within bright GRBs which in- 
tends to complement the work mentioned above, by specifically 
aiming at: 1) compare the results of the analysis of the spec- 
tra averaged over major pulses in the burst lightcurve with the 
time resolved spectra of the very same burst in order to quan- 
tify systematic differences; 2) consider both empirical (Band 
model) and more 'physical' (synchrotron shock model) spec- 
tral models for all spectra and compare the quality of the cor- 
responding fits. In particular each spectrum is fitted with the 
four models we choose (Band's, broken power-law, thermal 
Comptonization and synchrotron shock model, the latter fitted 
to temporal resolved BATSE spectra for the first time). We also 
examine any spectral 'violation' (with respect to the predicted 
slope in the case of synchrotron emission) for the entire burst 
evolution. The development of this work is the analysis and in- 
terpretation of the spectral evolution morphologies of this sam- 
ple of bursts. 

The paper is structured as follows: In sect. 2 we describe 
the data and their selection criteria, while the spectral mod- 
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els adopted for the analysis are detailed in sect. 3. Section 4 
presents the results of our work for the time integrated and 
time resolved spectra. Conclusions are drawn in sect. 5. The 
time evolution of the spectral parameter will be the content pre- 
sented in a following paper (Ghirlanda et al. in preparation). 

2. Data selection and analysis 

2. 1. Instrumental summary 

The Burst and Transient Source Experiment (BATSE) con- 
sisted of eight detection modules on the corners of the Compton 
Gamma Ray Observatory (CGRO - deorbited in summer 
2001). Each module was composed by two instruments: the 
Large Area Detector (LAD) designed for burst location and 
temporal analysis, and the Spectroscopic Detector (SD) suited 
for spectral study. 

Each LAD consisted of a thin circular Nal layer of 2025 
cm 2 collecting area. The nominal energy coverage was 28 - 
1800 keV (with minor variations between different detector 
modules). The energy resolution was about 20% at 511 keV 
(Fishman et al. 



1989a). 



The SDs were smaller (127 cm 2 ) but thicker and thus had 
a greater energy conversion efficiency. Their sensitive energy 
range varied with the gain settings and energy thresholds of 
the photomultipliers, but typically extended from 20 keV to 2 
MeV (Band et al. |1992[). The energy resolution is 7% at 662 



keV (Fishman et al. 1989a, Band et al. 1993) 



After a burst trigger, during the following 4 minutes burst 
mode, the LADs accumulated, among other data products, the 
High Energy Resolution Burst (HERB) spectra. These data are 
a sequence of 128 quasi-logarithmically spaced energy chan- 
nel spectra with a maximum time resolution of 128 ms. They 
were accumulated from the 4 most illuminated detectors. The 
SDs produced similar spectra with 256 energy channels and 
128 ms resolution (SHERB). The maximum number of spectra 
accumulated is 128 and 192 for the HERB and SHERB data, 
respectively. The LAD and SD spectral time integration algo- 
rithm was based on a count rate criterion so that the instruments 
provided spectra with a minimum integration time of 128 ms 
only for particularly intense bursts or around the peaks, while 
in most cases the spectral accumulation timescale was much 
greater than 128 ms. 

2.2. Data Selection 

The HERB data have been systematically preferred for this 
work because the higher detection area of the LAD ensures an 
higher count rate than the SD detectors. Despite their moderate 
energy resolution (if compared to the SDs ) they are suited for 
the continuum spectral study (Preece et al. 2000 ). 

There are a few cases (see col. 3 in Table |lj) for which the 
SHERB data were used because there were telemetry gaps. 

We have analyzed the HERB data from rank 1 detector. The 
rank of the detector (col. 4 in Table [j]) is an indication of the 
relative count rate during the trigger: rank 1 detector has the 
highest count rate, the best S/N and the highest spectral time 
resolution. For the SHERB data the detector choice is a com- 



promise between the highest degree of illumination (the first 
rank SD detector) and the highest gain which depends on differ- 
ent ground setting parameters, which have been changed during 
the mission (Band D., private communication). The gain of the 
photomultiplier, in fact, scales the spectrum up or down in en- 
ergy: the highest is the gain the more the detector is sensible to 
the low energy photons and the lower is the low energy thresh- 
old (Preece et al., |1996j ; Kitchin [l99"l| ). We selected the highest 
rank SD detector with the 511 keV calibration line centroid 
above uncompressed channel 800. 

2.3. The bursts sample 

The bursts were selected from the BATSE 4B catalog^] which 
is complete until 29 Aug 1996 (trigger 5586). This catalog was 
complemented with the on line catalog which includes the trig- 
gers after 5586 until trigger 8121 (9 Sept 2000). 

We have selected GRBs with a peak flux on t he 64 
ms time scale (calcu lated according t o Fish man et al. 1994, 
Meegan et al. |1994 Paciesas et al. |1999| ) higher than 20 
phot cm~ 2 sec -1 . This choice is motivated by the fact that 
bright bursts should provide time resolved spectra with good 
S/N (also on integration of the order of 128 ms). We collected 
in this way 38 bursts whose trigger number, name and peak flux 
are reported in Table, [j] (cols. 1, 2 and 6, respectively). 

The sample selected according to the peak flux criterion 
was reduced during the spectral analysis because of data prob- 
lems (trigger 1609, 1711, 3480) or because the number of spec- 
tra available for the spectral analysis was less than 5 (trigger 
1997, 2151, 2431, 2611, 3412, 571 1, 5989, 6293, 6904, 7647). 
This happened particularly in short bursts which have typical 
duration lower than 1 second. The bursts which were not ana- 
lyzed have dashes entries in tab [j]. The final sample of analyzed 
bursts contains 25 bursts and is obviously not complete with 
respect to the flux selection criterion. All the analyzed bursts 
except trigger 5563,7301,7549 are also present in Preece et al. 



(200C) spectral sample. 



3. Spectral analysis 

The spectral analysis has been performed with the software 
SOAR v3.0 (Spectroscopic Oriented Analysis Routines) by 
Ford (1993). The power of this software is its multi-level rou- 
tine structure which can be easily modified by the user: we 
implemented the spectral library with the broken power law 
and added the possibility of fitting each time resolved spectrum 
with a different energy binning scheme. The main capabilities 
of SOAR are: the possibility of analyzing HERB and SHERB 
data and rebinning (in time) spectra in different fashions (ac- 
cording to a fixed number of spectra per group or with a S/N 
ratio criterion). It is also possible to simulate burst spectra and 
spectral evolution (Ford 1993). 



We used the package XSPEC, as widely used for high en- 
ergy data, to test, for two bursts, the spectral analysis proce- 
dure. We found, within parameters errors, similar results. Any 



http://cossc.gsfc.nasa.gov/batse/4Bcatalog/index.html. 
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Table 1. The GRB sample 



Trigger^ 


GRB 


Data 


Detector 


P eakFlux c 


Background 


Peak c 


Spectra 1 


S/N s 






Type h 


TV UTTl 




[phot/ cm /sec] 


7r 


n d 






> 


143 


910503 


HERB 


LADO 


1 


52.08 ± 1.43 


8 


3 


4 


49 


— 


1473 


920311 


HERB 


LAD5 


1 


25.31 ± 0.7 





2 


3 


55 




1541 


920406 


SHERB 


SD2 


2 


38.32 ± 0.89 


5 


2 


1 


9 




1609 


920517 








67.59 ± 1.21 












1625 


920525 


HERB 


LAD4 


1 


28.06 ± 0.74 


8 


4 


3 


25 


— 


1711 


920720 








21.7 ± 0.7 












1997 


921022 








40.33 ± 0.85 






— 


— 


— 


2083 


921207 


HERB 


LADO 


1 


46.55 ± 0.92 





2 


2 


31 


— 


2151 


930131 








167.84 ± 2.63 






— 


— 


— 


2329 


930506 


SHERB 


SD2 


4 


42.57 ± 0.90 


11 


4 


1 


42 


— 


2431 


930706 








43.83 ± 0.89 






— 


— 


— 


2537 


930922 


HERB 


LAD1 


1 


27.28 ± 0.7 


31 


4 


3 


23 


45 


2611 


931031 








35.05 ± 0.89 












2798 


940206 


HERB 


LAD3 


1 


24.19 ± 0.73 


7 


3 


1 


16 




2831 


940217 


HERB 


LADO 


1 


44.34 ± 1.08 


10 


4 


3 


28 


— 


3412 


950211 








54.82 ± 0.76 






— 


— 


— 


3480 


950325 








21.61 ± 0.51 








— 


— 


3481 


950325 


HERB 


LAD2 


1 


25.7 ± 0.6 


6 


4 


3 


27 




3491 


950403 


HERB 


LAD3 


1 


30.65 ± 0.62 


6 


3 


3 


44 


45 


3492 


950403 


HERB 


LAD5 


1 


61.44 ± 0.91 


6 


4 


1 


22 


45 


3523 


950425 


HERB 


LADS 


1 


21.81 ± 0.53 


4 


2 


2 


50 




4368 


960114 


SHERB 


SDO 


1 


58.61 ± 0.83 


8 


4 


3 


29 


45 


5477 


960529 


HERB 


LAD1 


1 


29.35 ± 0.70 


7 


4 


2 


12 


— 


5563 


960804 


HERB 


LADA 


1 


22.7 ± 0.59 


7 


4 


1 


8 




5567 


960807 


HERB 


LADO 


1 


22.8 ± 0.6 


6 


4 


5 


25 


45 


5614 


960924 


SHERB 


SDO 


1 


183.37 ± 1.62 


12 


4 


1 


22 




5621 


961001 


HERB 


LAD2 


1 


25.64 ± 0.61 


31 


4 


2 


17 




5704 


961202 


HERB 


LADO 


1 


43.93 ±0.76 


26 


4 


1 


13 




5711 


961212 








41.25 ± 0.71 












5989 


970201 








77.61 ± 0.89 












6198 


970420 


HERB 


LAD4 


1 


66.5 ± 0.9 


12 


4 


1 


38 




6293 


970704 








88.53 ± 1. 












6404 


970930 


HERB 


LADO 


1 


24.01 ± 0.52 


21 


4 


2 


22 




6581 


980125 


HERB 


LADO 


1 


40.91 ± 0.69 


15 


4 


1 


17 




6904 


980706 








54.62 ± 0.87 












7301 


990104 


HERB 


LAD7 


1 


86.53 ± 0.93 


12 


4 


5 


24 




7549 


990506 


HERB 


LAD7 


1 


25.12 ± 0.58 


f) 


4 


6 


52 




7647 


990712 








24.06 ± 0.51 













a Burst reference number from the Gamma Burst Catalog at 

http : 1/ gammaray. msfc.nasa.gov /batse/ grb/ catalog / 'current / 

b (S)HERB: (Spectroscopic) High Energy Resolution Burst data. 

c Peak flux on the 64 ms time-scale, integrated over the 50 - 300 keV energy range. 

d Background is calculated on " # " number of spectra and with a polynomial function of degree "n". 

° Number of pulses spectroscopically analyzed within each burst. 

f Number of time resolved spectra per burst. 

g S/N method used in grouping the time resolved spectra accumulated by the instrument. 



difference can be ascribed to a different background calculation 
technique. 

The spectral analysis procedure consisted in the inspection 
of the light curve, i.e. the count rate in cts s _1 over the nominal 
energy band (typically 28-1800 keV) with each time bin corre- 
sponding to a time resolved spectrum (with typical integration 
time of 128 ms). 

The background is calculated on a time interval as close 
as possible to the burst, but not contaminated by the burst it- 
self. Typically we use a 1000 sec interval before and after the 
GRB, and fit the spectra contained in this interval with a poly- 
nomial background model examining the residuals for any sys- 
tematic problem. The background spectrum to be subtracted is 
obtained by extrapolating the polynomial best fit spectrum, for 
every channel, in the time interval occupied by the burst. As 
reported in Table [j], in most cases the degree of the polynomial 
is 4. 



The GRB light curve was used to identify the pulses as 
those structures in which the flux rises above and decays back 
to the background level. We accumulated the time resolved 
spectra within this time interval: this defines the time average 
pulse spectrum. In those pulses composed by several substruc- 
tures we treated them separately only if their flux starts and 
returns to the background level and if they are separated by 
at least a 128ms time bin. In the example trigger of Fig. [I] 
we accumulated two average pulse spectra, corresponding to 
the dashed regions in the plot. The time averaged pulse spec- 
trum was then fitted with the 4 spectral models described in 
the following section. The best fit results of the average spec- 
trum were used as initial parameters guess for fitting the time 
resolved spectra. 

For the time resolved spectra we used, as integration time, 
at first the shortest time available (i.e. 128 ms, at best), again 
limiting our analysis to the same time-span selected for the 
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Fig. 1. Example of light curve (trigger 2083). Each time bin 
represents the count rate summed over all the detectors that 
triggered, integrated in energy from 28 to 1800 keV and over 
the detection area (for details see Preece et al. 2000). The line- 



filled regions show the time intervals of the average pulse spec- 
trum accumulation. The start and stop times were chosen to 
have the same flux level. Our analysis is limited to the dashed 
areas, and therefore excludes the long slowly decaying part of 
the burst. 

time integrated analysis. The best fit parameters were then ex- 
amined in search for any indetermination: when at least 3 sub- 
sequent spectra had one or more undetermined spectral param- 
eter, the spectra were grouped according t o a S/N criterion. We 
followed the prescriptions of Preece et al. ( 200C ) who accumu- 
late subsequent HERB spectra until the S/N (calculated for the 
28-1800 keV energy range) is greater than 45; for the SHERB 
data the S/N threshold was fixed at 15 as already done by Ford 
et al. (1995). If after this accumulation the number of time re- 



solved spectra was less than 5 the burst was not included in the 
final list of analyzed events (see sect. 2.2). 

The spectra were also rebinned in energy in order to be con- 
fident that the Poisson statistic is represented by a normal dis- 
tribution in every channel for the application of the \ 2 mini- 
mization technique. We fixed the minimum number of counts 
per energy bin at 30 and 15 for the HERB and SHERB data, 
respectively. 

3.1. Spectral Models 

We used four spectral representations to fit the GRB time av- 
eraged and time resolved spectra. These functions were chosen 
in order to find whether a specific model can be considered as 
best representation of the spectral characteristics of pulses in 
bright GRBs. We explicitly exclude a fifth spectral model pro- 
posed in the literature by Ryde (1999) because it has a higher 
number of free parameters. 

3.1.1. The Band model 

This empirical model (BAND hereafter) was first proposed 
by Band et al. (1993) and, as already mentioned, is a good 



general and simple description of the time averaged (Band 
et al 



1993) and the time resolved spectra (Ford et al 



1995 



Preece et al. 1998). It contains the two continuum components 
in the keV-MeV band already discovered before BATSE: (i) 
a low energy power law with an exponential cutoff N(E) oc 
E a exp(—E/Eo) and (ii) a high energy power law N(E) oc 
E? (Matz et al. 119851). In fact the BAND model consists of 



2 power laws joined smoothly by an exponential roll-over, 
namely: 



N(E) 



A 



100 keV/ 



AEf> 



{a-l3)E Q 
100 keV 



exp 

a-f3 



for E < (a - j3) E 
exp(/3 — a) 

for E>(a-(3)E 



(1) 



where N(E) is in photons cm 2 s 1 keV 1 . The free parame- 
ters, which are the result of the fits, are: 

- A: the normalization constant @ 100 keV; 

- a: the low energy power law spectral index; 

- (3: the high energy power law spectral index; 

- Eq: the break energy, which represents the e-folding pa- 
rameter. 

If f3 < —2 the peak energy in the E Fe diagram (Fe is the flux 
in keV/cm _2 sec _1 keV) is E pea k = (a + 2)Eq and represents 
the energy at which most of the power is emitted. 

For the fitting procedure we had to assume an interval for 
the parameters and we fixed [—5,1] for a and (3 and the break 
energy was allowed to vary in the range [28-1800] keV. 

3.1 .2. The Broken Power Law model 

This model, called BPLW, is the simplest model used for fitting 
GRB spectra: it consists of two power laws sharply connected, 
with no curvature. Its analytical form is: 



N(E) 



AE a 

AE^E? 



for E < (a - 0) E 
for E>(a-f3)E 



(2) 



The free parameters are as before. In this model the peak energy 
of the EFe diagram coincides with the break energy Eq for 
P<-2. 



3.1 .3. The Comptonization model 

This spectral representation (COMP hereafter) is composed of 
a power law ending-up in an exponential cutoff, thus fitting 
well those spectra with a very steep high energy decline: 



N(E) = AE a exp ( - — 

Eq 



(3) 



The free parameters are as in the previous models without the 
high energy component. In fact the COMP model can be ana- 
lytically obtained from the BAND model assuming = — oo. 
This model might be considered to mimic the spectrum result- 
ing from multiple Compton emission from a thermal medium. 
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Fig. 2. Different electron energy distributions, solid line=post 
shock Maxwellian with a high energy power law tail (<5 = 3.5); 
dotted Z/«e=pre-shock relativistic Maxwellian. The distribu- 
tions are functions of 7/7*, and AT (7/7*) is in arbitrary units. 



3.1 .4. The Synchrotron Shock Model 

The model (SSM) is based on optically thin synchrotron emis- 
sion from relativistic p article s (either electrons and/or electron- 
positron pairs) (Tavani 1996 ). The electron energy distribution 
N(-f), which is a relativistic Maxwellian before the shock oc- 
curs, is modeled for the internal shock passage by adding a high 
energy powerlaw component (Tavani 1996): 



TV (7) oc 



3 

7" 



exp 



(-*) 



for 3- < 1 

7* — 

forl<i<2= 
— 7* — 7* 



(4) 



where 7m e c 2 is the electron energy, 7* = ksT* /m e c 2 is 
the pre-shock equilibrium electron energy, r y~ s is the supra- 
thermal component and j max is the maximum electron energy. 
The power law part of Nfa) (of index S) is related to the high 
energy spectral power law component of the photon spectrum 
(oc E p ) by the relation (3 = -{5 + l)/2. 

In Fig. H we report two cases of electron energy distribu- 
tion. Notice that in the case described by Eq. Q (solid line), the 
distribution extends on a wider energy range compared with 
the relativistic Maxwellian. In the figure a typical intermediate 
value for the power law index 5 = 3.5 is assumed. 

The synchrotron spectrum emitted by such an electron pop- 
ulation c an be obtained via the classic formulse(Rybicki & 
Lightman[i979|): 



I syn oc B ps ( / x 2 exp{-x) F 5/3 (-^) dx + 
Jo x 



e .11 



x ~ 5 F s/3( Jr) d.r 1 



where x = 7/7*, y = E/E* h and E* h is the synchrotron 
characteristic energy corresponding to those electrons with en- 
ergy E+ = rrieC 2 ^!*. B ps is the post-shock magnetic field 
which is the pre-e xistin g B enhanced by the strong shock pas- 
sage (Kennel et al |1984|)- 



Fig. 3. Synchrotron spectra corresponding to the two elec- 
tron energy distributions of Fig. ||: solid line=post shock 
Maxwellian with a high energy power law tail (S = 3.5); dot- 
ted //ne=pre-shock relativistic Maxwellian. The distributions 
are functions of vjv-k. The shown spectra are normalized to 
their respective peak flux. 



The shape of the synchrotron spectrum emitted by this elec- 
tron population (Fig. Q, solid line) has a substantial contin- 
uous curvature and is characterized by a low energy compo- 
nent (namely for E < E pea k ~ 2.5S*, where E+ = m^c 2 ^* 
is the electron energy) which should be well represented by 
the single electron synchrotron spectrum dependence ~ E 4 / 3 
(in the EFe diagram, i.e. ~ E^ 2 / 3, in the count spectrum) 
and by £ ,_ ( (5_3 )/ 2 at high energies (i.e. above E pea k). It is 
evident that this spectral form can account for the high en- 
ergy spectral variety of GRB observed spectra but has a fixed 
(—2/3) slope at low energies. One immediate prediction of this 
model is the relation between the electron energy 7*, the post- 
shock magnetic field B ps and the photon spectrum peak energy 
E* h — (3 e h/4ir m e ) B ps ^ 2 , which implies that the product 
< Bps'y 2 > can be constrained by the data. 

The free parameters of the model are the high energy elec- 
tron spectral index 6 and the characteristic energy E+ which 
describe the electron energy distribution. In the fit of this model 
to the time resolved spectra we assumed that the maximal in- 
terval of variation of these 2 parameters is 6 G [1, 10] and 
E* e [10 - 2000] keV. 



3.1 .5. Comparison of the fitting functions 

For clarity, in Fig. |] we report the spectral shapes correspond- 
ing to the four models considered. The potential of each model 
to better represent a particular spectral shape is clearly visible. 
The four models are similar at low energies but they differ in 
the high energy tail. In general we expect the average COMP 
model to overestimate the spectral break energy due to the lack 
of a power law high energy component. Furthermore the aver- 
age BAND spectrum would tend to be harder at low energies 
and softer at high energies than the BPLW model, due to the 
sharp break of the latter. The SSM average spectrum results 
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Fig. 4. The 4 spectral models: the spectral parameters corre- 
spond to the average values of the pulse time integrated spectra 
reported in Table ^. 

harder than the BAND and COMP model as the low energy 
spectral index is fixed in the SSM model, making the average 
spectral shape harder. 

We tested all the 4 spectral models presented above in or- 
der to be consisten t with the previous work present in literature 
(e.g. Preece et al. 2000 ). We notice, anyway, that the BPLW 
model is unphysical due to its sharp spectral break and the 
COMP model can be considered a subset of the more general 
BAND form, leaving the latter and the SSM as the primary 
spectral models to be compared. 

4. Results 

We have fitted the 25 GRBs of the sample with the four spec- 
tral models presented in sect. 3.1. For each peak present in the 
burst light curve we have analyzed the time integrated spectrum 
and the sequence of time resolved spectra. As mentioned in the 
Introduction the goals are the following: 

• determine whether there is a preferential shape which bet- 
ter represents the time integrated and/or the time resolved spec- 
tra of bright bursts; 

• alternatively, to determine the robustness of the spectral 
quantities with respect to the specific model considered; 

• quantify the difference between time integrated and time 
resolved spectra of the very same peak; 

• verify the strongest prediction of the synchrotron theory 
of a low en ergy s pectral shape not evolving and fixed to —2/3 
(Katz et al., Il994|)- 



For clarity we first present the results for the time integrated 
(sect. 4. 1), then for the time resolved (sect. 4.2) spectra, and fi- 
nally discuss the violation of the synchrotron 'limit' (sect. 4.3). 

4. 1. Time integrated pulse spectrum 

In Table ^ we report the best fit parameters of the 4 spectral 
models obtained by fitting the time averaged spectrum of the 
different peaks in each GRB. For multi-peaked bursts the num- 
ber of lines corresponds to the number of peaks analyzed. There 



are some gaps in each model columns corresponding to those 
pulses which are not fitted (i.e. the model parameters are unde- 
termined) by that model. In Fig. |]is reported an example of fit: 
the average pulse spectrum of trigger 3492 and the best fit and 
residuals for each model are displayed. 

4.1 .1 . Comparison of the spectral models 

From the time integrated spectral analysis we do not find a clear 
indication of a preferred fitting function to represent the spec- 
tra. In fact all the models give acceptable fits, and their xl e d 
are around one for all the 4 models although their median is 
definitely grater than one. However we can report of an indi- 
cation that the BAND model is the spectral shape which better 
fits the time integrated spectra. The BAND model has an av- 
erage Xred ~ 1-3' t0 be compared with the values 1.67, 1.63, 
1.74 of the BPLW, COMP and SSM0 models, respectively, and 
the width of the distributions, in terms of standard deviations is 
0.36 for the BAND and ~ 0.68 for the other three models. This 
indicates that the average pulse spectrum is better represented 
by the BAND model. We performed a maximum likelihood test 
and obtained that the improvement of the x 2 passing from the 
COMP model with 3 free parameters to the 4 free parameter 
BAND model corresponds, in most cases, to a better fit, al- 
though there exists cases where they are statistically indistin- 
guishable. This result agrees with what found by Band (1993). 
Many exceptions exists anyway: as an example let just mention 
the case of a spectrum whose high energy part is too steep to be 
represented by the BAND model power law (i.e. (3 = —5, first 
pulse of triggers 4368, 5567) while the COMP model can ac- 
commodate this fast decreasing spectrum with its exponential 
cutoff. Notice that in these cases also the BPLW model gives 
acceptable fits for (3: this is probably due to the sharp struc- 
ture of this model which makes the fit to overestimate the high 
energy spectral hardness (i.e. less negative (3 values). 

Note also that a statistically good x 2 e d °f course cannot be 
considered a sufficient condition for the quality of the fit and 
an analysis of the distribution of the residuals is necessary. 



4.1 .2. Spectral parameters distributions and 
average spectral shape 

We first computed for each model the distributions of the best 
fit parameters and found that they agree with t he res ults pre- 
sented in previous spectral studies (Band et al. 1993 ). In par- 
ticular for the BAND model we find that there is no corre- 
lation between the low and the high energy spectral indices 
(the Pearson's correlation coefficient is —0.006) but they group 
in the ranges -1.8 < a < -0.6 and -3.5 < (3 < -1.5 
which correspond to the intervals reported by Band et al. 
( 1993 ) and include the average val ues < a >— —1.03 and 
< (3 >= —3.31 found by Fenimore ( 1999 ). The peak energy of 
the BAND model E pea k ~ 205 ± 3 keV is instead harder than 



The SSM model has to be considered however apart. In fact due to 
the considerable number of average peak spectra which present an un- 
determined high energy supra-thermal index (8 > 10), the goodness 
of the model cannot be easily quantified. 
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100 1000 100 1000 

Energy (keV) Energy (keV) 

Fig. 5. Trigger 3492. Spectral fits to the pulse time averaged spectrum. The spectrum is integrated over the time interval [3.008- 
6.976] sec since the trigger time. The model fits and the residuals are for the BAND, BPLW, SSM and COMP model, displayed 
clockwise starting from the top left corner 



< E pea k >= 150±50 keV reported in Band et al. ( g993j )Q 
This is probably due to the fact that our burst sample repre- 
sents the bright end of the BATSE peak flux distribution, and 



The error associated with < E pea k > represents the spread of the 
values and the uncertainty on E pea k is the error on the average. 



because we integrated the time resolved spectra just around the 
peak, excluding the major decaying and rising parts of each 
pulse (see sect. ||) which are characterized by softer spectra. 

Table || reports the weighted average values of the best fit 
parameters of Table || for the pulse time integrated spectra. 
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Secondly, given that the analysis does not provide us with 
a best fitting function, let us consider the robustness of the pa- 
rameters found with respect to the choice of the spectral model. 

• The a distribution for the BAND model is peaked around 
— 1 and it is similar to the same parameter distribution for the 
COMP model; the BPLW, due to its sharp spectral break corre- 
sponding to the slope change E a — > E 13 , gives systematically 
lower a values and, in fact, its average is ~ —1.2. 

• The (3 distribution clusters around —2.3 in the BAND and 
BPLW model and for the latter the distribution is shifted to- 
wards higher values (i.e. harder spectra). The same happens 
with the SSM model having an harder average high energy 
spectral tail ((3 ~ -2.06) then the BAND model. 

• For the peak energy distribution we have that the weighted av- 
erage value for the BAND model is E pea k ~ 205 ± 3 keV. The 
BPLW break energy (= peak energy) is ~169 keV, whereas 
the COMP model, due to the lack of the high energy power law 
component, overestimates the spectral break, having a broad 
distribution and a weighted average of 233 keV. The SSM 
model peak energy is comparable with that of the BPLW model 
but lower than the BAND and COMP models. 

We conclude that the average spectral shape of the GRBs 
present in our sample does depend on the fitting model and the 
BAND and COMP model tend to give, although the latter lacks 
the high energy power law, comparable average spectral shape 
at low energies. 

4.1 .3. The Synchrotron limit violation 

The analysis of the low energy spectral indices shows that there 
are 4 pulses (those in italics in Table |J) whose average BAND 
best fit spectrum is harder, at 2a a level, than the limit E~ 2 ^ 3 
predicted by the synchrotron model. Three of these spectra, 
when fitted with the COMP model, show no a limit violation 
but their low energy power law spectral index (—0.667 ± 0.01, 
-0.64 ± 0.02 and -0.53 ± 0.08) is very close to -2/3. The 
same spectra have a poor fit with the SSM model showing sys- 
tematic trend in the residuals of the fits on most of the energy 
range. The average pulse spectra give only a weak indication 
of the violation of the synchrotron model low energy spectral 
index. On the other hand, as will be shown in the following, 
stronger indications of a violation come from the time resolved 
spectra. This occurs because the averaging of the spectral evo- 
lution over the duration of the pulse systematically yields softer 
spectra. 



Table 3. Weighted averages of the averaged-peak spectral pa- 
rameters for the four models. The break energy Ebrcak and the 
peak energy E poa k are in keV. 





BAND 


BPLW 


COMP 


SSM" 


<a> 


-0.977±0.003 


-1.185±0.002 


-1.013±0.002 


-2/3(/ix) 


</3> 


-2.27±0.01 


-2.09±0.004 




-2.06±0.06 


<E b reak> 


164.4±1.7 


161±0.7 


210.4±1.3 


67.4±0.7 


<E pca k> 


205±3 


161±0.7 


233±2 


168±2 



a For this model the spectral parameters are described in the text. 
Notice that f3 represents the high energy power law spectral index 
derived from the slope S of the particle distribution, as /3 — — (5 + 
l)/2. 




3 4 5 6 7 

time [sec] 



Fig. 6. Trigger 3492. Spectral evolution of the COMP model 
parameters fitted to the time resolved spectra. Light curve on 
the 64ms time-scale (panel a), for the energy range 110-320 
keV (corresponding to channel 3 of the detector); low en- 
ergy spectral index (b), dotted line: the average a value from 
Table]|, dashed line: the synchrotron model limit; peak energy 
(c), dotted line: the average Eb rea k of the COMP model. 



4.2. Time resolved spectra 

Let us now consider the results of the fitting of the time re- 
solved spectra with the same four models. 

In Fig. |6] an example of spectral evolution is reported: for 
each peak we obtain a sequence of best fit parameters relative to 
the time resolved spectra which characterize the temporal evo- 
lution of the spectrum during the burst. While a full description 
and discussion on the parameter evolution will be presented 
elsewhere (Ghirlanda et al, in preparation), we note that both 
the peak energy and the low energy spectral index evolve in 



phase: the flux and the spectrum becomes harder during the rise 
phase and s ofter d uring the decay, a behavior already found by 
Ford et al. ( 1995 ) as a characteristic spectral evolution mor- 
phology of several bursts. 

In the following we report and comment on the parameter 
distributions of the 4 models for a comparison (i) among the 
model themselves and (ii) with the results on the pulse average 
spectrum. Finally we examine the SSM limit violation in some 
well defined cases. 
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Table 2. The peak average spectral results: for each trigger (col. 1) the best fit parameters, of the four models described in the 
text, are reported: for multi peaked bursts all the peaks have their set of spectral parameters. The pulses which violate the low 
energy synchrotron limit are reported in Italics. Energy is in keV. 
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7549 










-1.19±0.01 


-1.76±0.()2 


247±11 


1 


.5 


-1.13±0.01 


1192±59 


1 


.4 


182±8 


2.6±0.05 


1.7 




-0.98±0.03 


-2.5±0.3 


494±46 


1.2 


-1.19±0.01 


-2.17±0.()7 


282±15 


1 


.4 


-1.01±0.02 


558±36 


1 


.2 


276±22 


4.0±0.3 


1.3 




-0.68±0.02 


-2.64±0.18 


397±23 


1.1 


-0.94±0.01 


-2.07±0.04 


262±8.3 


1 


.9 


-0.72±0.02 


446±18 


1 


.2 


650±82 


8.7±3.2 


1.7 




-0.97±0.03 


-2.2±0.07 


277±22 


1.2 


-1.28±0.01 


-2.08±0.04 


180±7 


1 


.8 


-1.09±0.02 


396±21 


1 


.6 


161±9 


3.8±0.2 


1.1 




-0.84±0.02 


-2.67±0.09 


201±8 


2.0 


-1.22±0.01 


-2.24±0.03 


156±3.1 


4 


.3 


-0.91±0.02 


232±7 


2 


.3 


231±15 


8.4±1.6 


1.7 



a Energy in keV. 



4.2.1 . Comparison of the spectral models 

For a general comparison on the quality of the fits with the 
different models we have plotted in Fig. || their reduced x 2 dis- 
tributions. These are all centered around one and again it is not 
possible to identify any preferable spectral model even con- 
sidering their spread. Anyway we can note some differences: 
considering that all these distributions are asymmetric towards 
2, we fitted a 6 parameter function, namely a right asymmetric 



gaussian distribution, and obtain that the BAND and COMP 
model have the lower dispersed distributions with a ~ 0.14, 
to be compared with — 0.2 for the BPLW and SSM model. 
This result indicates that in terms of reduced x 2 the BAND and 
COMP model could better represent the time resolved spectra 
of bright bursts. Also in this case some counter examples ex- 
ist showing that in general within a single pulse time resolved 
spectra can be fitted by different spectral models (see Fig. R). 
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We also tested, with a Kolmogorov-Smirnov test, whether the 
BAND and COMP distributions could have been drawn from 
the same distribution and obtain that this is the case with a prob- 
ability of 0.99 (95% confidence level). 



4.2.2. Time resolved vs time integrated spectra 



100 - 



Let us now compare with the corresponding results for the time 
integrated spectra. As an example in Fig. [7] we report the peak 
spectrum of the trigger 2083. Comparing with Fig || (showing 
the spectrum time integrated over the whole peak, for trigger 
4392) it is evident that the time resolved spectra better con- 
straint the best fitting model: in fact (for example in this case) 
the time averaged pulse spectrum was satisfactorily fitted by all 
the 4 models whereas the peak time resolved spectrum is not 
well fitted by the SSM model because the low energy power- 
law spectral index is fixed at -2/3 and the spectrum is harder 
than this slope (a = -0.03 ± 0.06 for the BAND model). 
In this example the BAND, COMP and BPLW model provide 
the best fits and among them the BAND model has the best 



= 0.99. 



This represents a clear indication that the time resolved 
spectra have to be used in determining the spectral properties of 
GRBs: average spectra can be effective in comparing the global 
properties of different bursts, but the actual spectral shape re- 
quires high temporal resolution data. 
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4.2.3. Spectral parameters distributions 

In this Section we consider the distribution of the model param- 
eters as inferred from the time resolved data, compare again the 
different models, and compare these results with the time inte- 
grated ones and with previous findings. 



Fig. 8. Reduced x 2 distributions for the 4 spectral models. The 
total number of time resolved spectral fits is 700 for the BAND 
and BPLW (dotted and solid line top panel) and COMP and 
SSM (solid and dot-dashed line bottom panel). 



The low energy spectral component. 
The BAND, COMP and BPLW fits are comparable at low 
energies as shown by the corresponding a distributions pre- 
sented in fig. ||: the BAND and COMP model distributions 
are similar, as in the case of the time integrated spectra, and 
both have a mode of —0.85 ± 0.1, consistent at 2er with the 
BPLW average value —1.15 ± 0.1. Note that qualitatively 
the extension of the a distribution of the BPLW model to- 
wards lower values could be attributed to the fact that at 
low energies the sharp break tends to underestimate the 
hardness of the spectrum compared to a smoothly curved 
model. 

The average low energy spectral slope obtained from the 
time resolved spectra is harder than what obtained with 
the time integrated pulse spectra for all the three models 
(BAND, BPLW, COMP). This is a consequence of time in- 
tegration (i.e. hardness averaging) of the spectral evolution 
(which can be also very dramatic) over the entire rise and/or 
decay phase of the pulse. 

Even though we present the spectral analysis with 4 mod- 
els separately we can co mpare these results with those ob- 
tained by Preece et al. ( 2000 ), where they considered the 
low energy spectral index distribution regardless of the 
(best) fitted model. We obtain a similar distribution but with 



a harder average low energy spectrum (a ~ —0.85 ± 0.1 
considering the BAND and COMP model, to be compared 
with ~ —1.05 ± 0.1). This difference might be a conse- 
quence of the fact that we are considering a subsample of 



that of Preece et al. (2000), including only the bright GRBs 



which might be intrinsically characterized by a harder spec- 



trum (Borgonovo & Ryde 2001 ) and because we restricted 



our time resolved spectral analysis to the peak interval ex- 
cluding the inter-pulse phase of multi-peaked events. 

The high energy spectral component. 
In many spectra the steepness of the count distribution 
above the break energy Eq, coupled with the poor S/N ra- 
tio, result in a poorly determined value of /3. This is true 
especially in the case of the BAND model: when the break 
energy E$ sets at high values the exponential roll-over 
leaves too few high energy spectral channels free for fit- 
ting the E 13 component properly. In Fig. we report the 
BPLW, BAND and SSM (3 distributions. Also in this case 
the BPLW model, due to its sharpness, tends to overesti- 
mate the hardness of the count spectrum giving systemati- 
cally higher values of (3 than the BAND model. The mode 
of this parameter is -2.45 ± 0.1 and -2.05 ± 0.1 for the 
BAND and BPLW model, respectively. 
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Fig. 7. Trigger 2083. Spectral fits to the pulse peak spectrum. The spectrum is integrated over the time interval [1.088-1.216] sec 
since the trigger time. The model fits and the residuals are for the BAND, BPLW, SSM and COMP model displayed clockwise 
starting from the top left corner. 



The SSM average (3 is —2.17 which is consistent with what 
found from the average pulse spectral analysis. The average 
f3 value for the BAND model (-2.45±0.1)is softer than the 
time averaged results, and the BPLW model value (—2.05) 
is consistent with the results of Table. [}J 



In fig. it is also indicated the critical value —2: spectra 
with j3 > —2 are rising in EFe, and we can set only a lower 
limit to Ep e ak- These cases correspond to the 18% and 7% 
of the spectra fitted with the BAND and BPLW model and 
the 16% for the SSM. Notice that there is also a subclass of 
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Fig. 9. Low energy power law spectral index (a) distributions 
derived from the time resolved spectral analysis. Solid line: 
BPLW model, dotted line: BAND model, dashed line: COMP 
model. The vertical line represents the synchrotron limit (a = 
—2/3) for the low energy spectral shape. 



Fig. 11. Peak energy distribution for the 4 spectral models. 
solid line: BPLW model, dotted line: BAND model, dashed 
line: COMP model, dot-dashed line: SSM model. Spectra with 
undetermined peak energy (i.e. the high energy threshold 1800 
keV assumed as lower limit) are reported in the last bin. 




-5 -A -3 -2 -1 

(BPLW, GRB, SSM model) 

Fig. 10. High energy power law spectral index {(3) distributions 
for the time resolved spectra. Solid line: BPLW model; dotted 
line: BAND model; dot-dashed line: SSM model, in this case 
(3 is calculated from the S parameter (see text). Also shown 
(bin with /3=-5) the time resolved spectra with undetermined 
high energy spectral index for the BAND model. The dashed 
vertical line indicates the critical value [3 = —2. 



In Fig.|ll| the peak energy distributions for the various 
models are reported and it is also shown the bin with 
E pea k=1800 keV, assumed as lower limit of the peak en- 
ergy for those spectra in which the BAND, BPLW or SSM 
model have (3 > — 2. 

The mode is E peak = 2801^ keV for the BAND model, 
consistent, within its error bar, with the BPLW most proba- 
ble value of 211±f| keV. The COMP model, instead, gives 
a highly asymmetric peak energy distribution with a mode 
of 595^8g 4 keV because the lack of a high energy power 
law component tends to overestimate the energy corre- 
sponding to the start of the exponential cutoff. The SSM 
model has an average E pea k ~ 316lg2 keV with a wide 
distribution. 

From the analysis of the average spectral shape, obtained 
from the time integrated and from the time resolved spec- 
tral analysis, we would like to stress that if we want to char- 
acterize the spectral hardness it is necessary to extend the 
single parameter analysis, typically based on the peak en- 
ergy or the hardness ratio, and consider the low and high 
energy spectral components. 



soft-spectra with j3 < — 3 which are characterized by a very 
steep spectral tail: these spectra are clearly better fitted by 
the COMP model. 

The spread in these (3 distributions (—4 < j3 < — 1.5) cor- 



responds to what found by Band et al. (1993) 



The spectral break. 

The most important spectral parameter obtained in fitting 
the spectrum with these models is E pea k- As just mentioned 
this characteristic energy can be obviously calculated only 
for those spectra (BPLW and BAND model) with [3 < -2. 



4.3. The Synchrotron limit violation 

A well known prediction of the optically thin synchrotron 
model is that the asymptotic low energy photon slope a should 
be lower than or equal to —2/3 (Katz et al. 1994). From the 
analysis of the time resolved spectra we obtain that not only 
the low energy power law slope can violate this limit, but it can 
also evolve dra matical ly during a single pulse (as already found 
by Crider et al. |l997a[ b). 

We can characterize this behaviour, for example, via the 
BAND and COMP model fits (we exclude the BPLW model 
which, as showed above, tends to underestimate the hardness 
of the spectrum at low energies). We obtain that the 14% of the 
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Fig. 12. Trigger 2083. Spectral evolution of the BAND model 
parameters fitted to the time resolved spectra. Light curve on 
the 64ms time-scale (panel a) for the energy range 110-320 
keV (corresponding to channel 3 of the detector); low energy 
spectral index and (dashed) synchrotron shock model limit 
a = —2/3 (b); peak energy (c). For reference the average val- 
ues of a and E pea k obtained from the time resolved spectra 
(dotted line) and the synchrotron model limit (dashed line) are 
reported. 

time resolved spectra fitted with the BAND model are incon- 
sistent with a < —2/3 at 2cr. A similar percentage (~1 1.7% of 
course mostly for the same spectra) of spectra violating the a 
limit is found for the COMP model. 

For these extremely hard spectra no correlation between a 
and any other fit parameter is found, and the violation occurs 
both in the rise and decay phase of the pulses. Moreover ~ 21% 
of the time resolved spectra corresponding to the peak time bin 
violate at 2a the synchrotron limit, indicating that this violation 
happens during the peak phase and not preferentially before or 
after it. As an example in Fig. 12 we show the spectral evolu- 
tion of a in the case of trigger GRB921207 (BAND model). 
During the main peak (t£ [0.0,4.5] sec) the majority of the time 
resolved spectra violates the synchrotron model limit (dashed 
line in panel (c)). Notice that during this time interval the low 
energy spectral index evolves between —2/3 and 0.0 and the 
peak energy (panel (b)) decreases monotonically. Another in- 
teresting example is reported in Fig. |l3[ GRB960924 shows a 
low energy spectral shape harder than —2/3 only around the 
peak, and the a evolution during the rise phase covers an inter- 
val of Aa ~ 1.4. 

The predictions of the synchrotro n mod el have been re- 
cently discussed by Lloyd & Petrosian fl200lj L&P2001). They 



Fig. 13. Trigger 5614. Spectral evolution of the COMP model 
fitted to the time resolved spectra. Light curve on the 64ms 
time-scale panel (a) for the energy range 110-320 keV (cor- 
responding to channel 3 of the detector); low energy spectral 
index and (dashed) synchrotron shock model limit a = —2/3 
(b); peak energy (c). For reference are reported the average val- 
ues of a and E pea k obtained from the time resolved spectra 
(dotted line) and the synchrotron model limit (dashed line). 

claim that the low energy spectral limit —2/3 is not so con- 
straining if (among other assumptions) an anisotropic pitch an- 
gle distribution is assumed for the emitting electrons (which 
should be the case in a low density and intense magnetic field 
regime). In the case of small pitch angles (\P <C 1) the emitted 
spectrum should be characterized by three components: a low 
energy fiat count spectrum t>° (for v <C v s = 2/3 Vb/ilm ^ 2 )) 
followed by the typical j/~ 2 / 3 power law (for v B < v < v m ) 
and then by ^~(p +1 )/ 2 . p is the electron power law energy in- 
dex, v s the break at the transition from the fiat spectrum to the 
—2/3 slope, and v m represents the frequency corresponding to 
7 m where the electron energy distribution is smoothly cut off 
(at low energies). 

According to L&P2001 a > -2/3 is allowed and the 
slope v a obtained from fitting a two component model (like 
the BAND function) to such a three components spectrum is 
typically an average between v° and v~ 2 ^ 3 . During the emis- 
sion the electrons cool so that their low energy limit j m de- 
creases and as a consequence the peak energy becomes softer 
Epeak oc 7^. The electrons average pitch angle decreases (al- 
though their distribution can be still anisotropic) and this causes 
the frequency v s oc ('j m 'I' 2 ) -1 to increase. The spectral evo- 
lution predicted by this model then is of hardening of the low 
energy powerlaw - because of the progressive disappearance 
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Fig. 14. Trigger 1625. Spectral evolution of the BAND model 
parameters fitted to the time resolved spectra. Light curve on 
the 64ms time-scale (panel a) for the energy range 110-320 
keV (corresponding to channel 3 of the detector); low energy 
spectral index and {dashed) synchrotron shock model limit a = 
—2/3 (b); peak energy (c). 



Fig. 15. Trigger 5477. Spectral evolution of the BAND model 
parameters fitted to the time resolved spectra. Light curve on 
the 64ms time-scale (panel a) for the energy range 110-320 
keV (corresponding to channel 3 of the detector); low energy 
spectral index and (dashed) synchrotron shock model limit a = 
—2/3 (b); peak energy (c). 



of the v~ 2 l 3 spectral component - while the peak energy natu- 
rally evolves from hard to soft. Thus in burst violating the syn- 
chrotron limit, a negative correlation between the peak energy 
and the low energy spectral index is expected. 

We report in Fig. [IJ as an example the evolution during 
GRB 920525 which shows a low energy spectral component 
harder than —2/3 during the main peak (for ~ 2 sec). We note 
(panel (b) and (c) in Fig. |l4|) that there is no evidence for neg- 
ative correlation between the peak energy and a: during the 
rise phase of the flux (panel (a)) the peak energy increases 
regardless of the rise and partial decay of the spectral index, 
while for the rest of the pulse the peak energy decreases and 
alpha varies above the —2/3 limit. The same happens for GRB 
921207 (Fig. |l§) and GRB 960529 (Fig. |l5j): in these bursts 
the peak energy decreases while a goes above and below the 
synchrotron limit. 

5. Conclusions 

We considered a sample of bright bursts detected by BATSE 
and performed a uniform analysis for the time integrated and 
the time resolved (typically 128 ms) spectra with four different 
models proposed in the literature. 

We find that even with this time resolution no parametric 
model can better represent the data and different spectra require 
different shapes, re-confirming the erratic behaviour of bursts 



and also possibly indicating that time resolution on time-scales 
comparable with the variability one is needed to shed light on 
such erratic characteristics. 

Indeed, an important result we confirm is that the aver- 
age time integrated spectrum often used in the literature does 
not well represent the very same event resolved on shorter 
time-scales. The time integrated spectra might still be used 
for a comparison of the average spectral shape among differ- 
ent pulses and as indicators of the average spectral parameters 
of the time resolved analysis although only the time resolved 
spectra should be used in any test of a physical emission model. 

Finally, a considerable number of the fitted spectra are char- 
acterized by extremely hard low energy components with spec- 
tral index a greater than —2/3, value predicted by synchrotron 
theory (Katz 1994). This violation was found by Crider et al. 
( 1 997a ) and has been recently reported by Frontera et al. ( 2000| ) 
in some GRBs observed by BeppoSAX. They report 1 sec time 
resolved spectra significantly harder than E~ 2 / 3 , mainly dur- 
ing the first phase of the burst emission. We have found that 
in 1 1 of the 25 bursts analyzed the a limit violation is evident 
mainly in the spectra around the peak both during the rise and 
decay phase, and this could indicate that at least at some stages 
of the burst evolution - possibly near the peak of emission itself 
- radiative processes, other than synchrotron, can dominate the 
emission. 
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We also reported some examples of bursts which violate 
the synchrotron limit and are not characterized by the a-E pea k 
anticorrelation predicted by the small pitch angle distribution 



synchrotron model proposed by Lloyd & Petrosian (2001 ). 

The obvious extension of this work, namely the study and 
discussion of the temporal evolution of the spectral shape, will 
be the subject of a forthcoming paper (Ghirlanda et al., in 
preparation). 
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